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Experiments have shown that the microwave ionization probability of a highly excited almost 
monodimensional hydrogen atom subjected to a microwave pulse sometimes grows in steps when 
the peak electric field of the pulse is increased. Classical pulsed simulations display the same 

^\ • steps, which have been traced to phase-space metamorphoses. Quantum numerical calculations 

again exhibit the same ionization steps. I show that the time-sequence of two level interactions, 
responsible for the observed steps in the quantum picture, is strictly related to the classical phase 

O^J ■ space structures generated by the above mentioned metamorphoses. 
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I. INTRODUCTION 



In classical physics, integrability - even for nondissipative systems - is the exception [1]. Regular - both periodic and 
quasiperiodic - and chaotic classes of motion can coexist in phase space. Classes of motion can appear or disappear 
when the system parameters are changed; periodic orbit bifurcation points and thresholds for chaos are among the 
most widely known of the critical points where this happens. 

The relevance of nonintegrable classical dynamics to quantum dynamics has been [2] and still is the source of much 
discussion, which has lately concentrated mainly on three problems: dynamical tunneling, both chaotic or resonance 
assisted [3] , the construction of atomic nondispersive wave packets [4] , and the ionization of Rydberg atoms [5] . 

Existing semiclassical theories predict a direct dependance of the quantum wavefunction on the corresponding 
classical physics, for instance via the classical action function (when it can be built). But the temporal regions of 
applicability of such theories generally are not well delineated, especially near the critical values and for chaos. Little 
is known about the case where a parameter is swept in time through a classical critical value. Under what conditions 
|' and to what extent can the evolution of a nonintegrable quantum system be "guided" by short-time nonadiabatic 
behavior present in its corresponding nonintegrable classical evolution? 

One semiclassical dynamical system which has been for a long time the object of extensive experimental study 
is a highly excited hydrogen atom exposed to a partially ionizing short pulse of microwave electric field. If the 
atom is prepared in a static electric field collincar with the microwave field, and an extreme Stark state is excited, 
a one-dimensional approximation can be used, thus greately simplifying the treatment. Experimental and quantum 
numerical results are often found to be nearly classical, agreeing reasonably well with numerical predictions based on 
classical models. The classical electron dynamics is a microwave frequency sensitive mixture of regular motion, locally 
\Q • chaotic motion, and globally chaotic motion [6]. The last motion eventually results in ionization of the classical atom. 
Resonant regular motion is displayed as island features in stroboscopic surfaces of section of the phase space. 

Let the initial principal quantum number of the atom be no- When the ratio ujq = n^oj of the microwave frequency 
u) to the initial Kepler orbit frequency l/rig is near unity, resonant classical motion dominates the character of a 
J> , number of semiclassical Floquet eigenstates which are experimentally populated at the peak of the pulse [6] . This 
inhibition of chaos produces an experimentally observed increase in the threshold field for ionization in the quantum 
system [7] . Bounding the island region in phase space is a locally chaotic region which can be crossed by the system 
during the rise and fall of the pulse [8] . Such crossings are examples of a phase-space transport process which would 
not occur in a fixed-field experiment [9, 10]. They provide an explanation for observed nearly classical final quantum 
state distributions which are double peaked as a function of the principal quantum number [11]. I shall call principal 
primary resonance zone the u>q = 1 resonance island together with the surrounding band of local chaos. 

For different constant values of the microwave field strength, the structure of the resonance zones can be different. 
The structural changes have been called phase-space metamorphoses [12], and it has been suggested that they could 
influence the microwave ionization process. For example, in Ref. [12] it was shown that the destruction of the last 
invariant curve between a primary resonant zone and the global chaos above it can lengthen the classical electron's 
escape time for ensembles of initial conditions in the global chaos region. Here I am interested in the metamorphoses 
of the internal structure of the resonance islands themselves. These metamorphoses are associated with periodic orbit 
bifurcations which generate, within the primary resonance zones, new "secondary" resonance Kolmogorov-Arnold- 
Moser (KAM) island chains, and with the destruction of these same secondary resonance island chains. 

Laboratory observations have been reported of sequences of steps in the (experimentally defined) ionization of 
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almost monodimcnsional Hydrogen Rydberg atoms versus the peak field strength of the microwave pulse [13] . Classical 
numerical simulations reproduce such steps. These classical step have been shown to be due to metamorphosis-induced 
variations of the probability for the transport in phase space responsible for ionization [13]. Together with more 
extensive numerical evidence and a more detailed analysis of the classical process [16], I now present an explanation 
of the good agreement of the quantum simulations with the classical ones: even if, in terms of Floquet states, the 
quantum evolution is dominated by a sequence of two level processes (avoided crossings traversed by the system) 
and can therefore be viewed as deeply quantum, the interactions between these couples of states exist because of the 
classical secondary island chains. 

The present paper is organized as follows: in section II I present the system and the numerical techniques I used to 
investigate it; the results of my investigation are presented in section III. Section IV is dedicated to the interpretation 
of these results; section V finally sums up my findings. 

II. THE SYSTEM AND NUMERICAL METHODS 

The system I investigate is a ID model for a stretched highly excited hydrogen atom in collincar static and 
monochromatic microwave electric fields [14]. In atomic units the Hamiltonian reads 

V 2 1 

H = — h z[F{t)sin(ujt + cf>o) — F s ], z > 0, (1) 

2 z 

where F(t) is the strength of the microwave field, 4>o is its initial phase, and F s is the static field strength. To simulate 
the experimental situation described in Rcf. [13], the envelope of the microwave pulse is chosen to be 

F(t) = F max sin(TTt/T p ), 

where T p = 7.78 ns is the length of the microwave pulse. For most of the simulations in the present paper this 
corresponds to 136 microwave periods. 

As is costumary for this system I shall use classically scaled parameters [14]: F™ ax — n\F max , F s0 = tiqF s , and 
uo = tiqUo. The scaled frequency is furthermore corrected to the first order in F s0 to compensate for the Stark shift 
of the atomic frequencies [6] : 

cj' = w /(l - 3F s0 ). 
A. Classical methods 

Following Rcf. [15], classical numerical integration in Ref. [16] was performed in free-atom action angle variables 
(/, 9), valid when the electron's energy is negative and defined by the equations: 

z = 2I 2 sin 2 (£/2), 
p=(l/J)cot(£/2), 

where the eccentric anomaly angle £(0) is defined by 9 = £ — sin£. 

To avoid equations of motion containing terms which diverge as z approaches zero, a dummy time 77, defined by 
dt = (1 — cos £)dr) was also introduced [17]. The true time t increases monotonically with r\. Newton's equations of 
motion then become 

dl/dr) = -I 2 sm£[F(t) sin (tot + <j> ) - F s ], 
d^/dr] = 1/I 3 + 2/(1 - cos£)[F(i) sin (ut + <fo) - F s ], 
dt/dr/ = 1 — cos£. 

The above three ordinary differential equations were numerically integrated using a fixed step, fourth order Rungc- 
Kutta routine. The ensemble of initial conditions (9o,Iq) was chosen so that, to first order in F s , the (classical) 
electron energy E(t — 0) be equal to the value of the energy of the experimentally prepared initial quantum energy 
eigenstate: 



~2^ - l F ° n ° = E ^ = ~27f ~ Fs/ ° 2(1 " COS€o) ' 



(2) 
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Given a value of #o(£o)j this equation determines a corresponding value of Iq. 

Classical values for the "ionization" probability Pj at the end of the pulse were averaged over uniform distributions 
of the initial angle 80 and microwave field phase 4>q. The "ionization" probability contains two contributions. One 
comes from trajectories which are terminated at some time during the pulse where the instantaneous value of the 
energy E exceeds the value — 2\/77 f° r rapid ionization in the static field alone. To match the experimental definition 
of "ionization" [13], the other contribution comes from trajectories whose final value of E corresponds, according to 
equation (2) above, to energy eigenstates with quantum numbers outside the interval n € [50, 90]. 

Stroboscopic surfaces of section of the long-time evolution in (0,1) space were also computed for F(t) equal to 
various constants F [16]. These reveal the long-time structures existing in phase space for the various values of F 
instantaneously traversed during the microwave pulse. For a constant amplitude microwave field, the action-angle 
generalized phase space (in the three coordinates /, 8, and t) is a well known [6] admixture of zones of regular and 
chaotic motion. Zones of regular electron motion include those occupied by nested vortex tubes and those of the 
field- modulated atom [8]. Chaotic zones include zones of local chaos and of global chaos. Stroboscopic surfaces of 
section (Poincare maps) reflect these zones in the full phase space as characteristic zones in the action-angle (J, 8) 
subspaces ("sectioned" phase space): the vortex tubes produce resonance island zones and the local chaos produces 
zones of irregular motion which have been called "separatrix" zones [18] (the actual separatrix of integrable systems 
is "broken" in the nonintegrable driven hydrogen atom problem and becomes a chaotic layer which grows with 
increasing microwave field strength). These structures are self-similar: secondary resonance zones are found within 
primary resonance zones and so on; but while all primary resonance zones are born at zero microwave field, periodic 
orbit bifurcations produce secondary resonance zones at nonzero field values. 

The comparison of instantaneous ensemble distributions in phase space with these long time surfaces of section, 
for various times t during the pulse [13] has proven a useful tool for the understanding of the evolution of the 
ensemble itself: when a phase space structure influences the pulsed evolution, the ensemble appears deformed by it in 
a characteristic way. For example, an ensemble of points initially on a straight line is stretched by a resonance island 
into a whorl [19] as the points of the ensemble move along the island's invariant curves. Local whorl-like deformations 
of the ensamble are therefore indication of the influence of an underlying island which can be identified with the help 
of the surface of section at that time. Similarly, local hyperbolic chaos deforms the ensemble into tendrils [19]. 



B. Quantum methods 



The time evolution of the quantum system is evaluated by numerical integration of the Schrodinger equation with 
H given by eq.(l) on a finite subset of the (bound) free atom basis ip n (z) chosen big enough so that the probability 
reflected by the borders is small [15]. Given an initial state tp(z,0) — S n C n (0))/i n (z), the equation for the evolution 
of the expansion coefficients C n (t) is 

i — — E n C n (t) + T(t)Y, m Z ntm C m (t), 

T{t) = (F(t)sin(cjt + <f> ) - Fs), (3) 

where Z n ^ m is the matrix element of the operator z between the states n and m. We approximate !F(t) with the 
function !F'(t) — J r (t)At'Ek^(t — kAt) which tends to T(t) for At — > (in the sense that the integral of their difference 
over an arbitrary time interval goes to zero as At) [15]. This procedure involves unphysical parameters which have 
to be carefully chosen as not to falsify the results of the integration: I discuss my choices in Appendix A. 

The initial state for our simulations is taken to be the eigenstate with quantum number no of the atomic Hamiltonian 
in the static field alone, calculated in the chosen finite subset of the free atom basis. 



1. Quasienergy curves 

Like in the classical case, an understanding of the dynamics at constant microwave amplitude will help illuminating 
the pulsed dynamics. If F(t) is a constant, Schrodingcr's equation is a differential equation with time-periodic 
coefficients. Floquet's theory is therefore applicable [20] and tells us that the equation has solutions in the form 

^ (t) = (t)e~^\ ^ (t + T) = ^ (t) 

where T is the period of the microwave, the constants Si take the name of quasienergies and the functions ipi(t) are 
called quasienergy (or Floquet) states. For periodic systems quasienergies take the place of energies in the description 
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of the dynamical properties of the system. It can be easily seen that the same physical state ipi(t) can be written as 

^(f) = \^.(^ e -i^nt/T^ e -ie it /h+i2 V nt/T 

where $>i(t)e~ l27Tnt / T is again a T-periodic function. We therefore have an infinity of replicas of the quasienergy 
spectrum shifted by 2irh/T one respect to another. It will then be sufficient to restrict ourselves to the energy interval 
[0, 2nH/T) (first Brillouin zone) to have all the levels. In particular, since the free atom has a continuum, this too will 
be brought into the first Brillouin zone. Due to the presence of the free atom continuum, the microwave transforms 
all the free atom states into resonances with finite widths I\: the quasienergy spectrum for an atom in microwaves is 
therefore absolutely continuous (no point spectrum exists) [21]. 

For the actual calculation of the quasienergies it is useful to resort to the time evolution operator G over one 
period T of the microwave (Floquet operator). There isn't just one Floquct operator, but a whole family of them, 
parametrized by the phase <j> = (uito + <f>o) of the field at the beginning of the period. While the cigenfunctions $i((j>) 
of each member of the family are different and represent the different spatial structures of the states at different 
times during the microwave period, the family shares the same eigenvalues Gi. The eigenvalues Gi of the Floquet 
operator and the quasienergies are connected by the relationship Gi — e ^ lSiT / h . One common procedure to obtain the 
quasienergies is therefore to numerically calculate a one period evolution operator and diagonalize it thus obtaining 
the GiS. The above relationship then gives us the quasienergies on the circle [0, 2irh/T). 

To investigate the evolution of a pulsed system, and also -as we shall see in Appendix A- to control numerical 
errors and approximations, it is useful to plot the quasienergies as functions of the microwave field strength F. The 
resulting quasienergy curves (see e.g. Fig. 1) help us connect the characteristics of the system at different values 
of the microwave field strength. The quasienergies change with F and undergo a sequence of avoided crossings of 
different widths. Since in the II? hydrogen atom all the quasienergies fall into the same symmetry class [22], the Von 
Neumann- Wigncr theorem [23] tells us that they all repel each other; this implies the absence of real level crossings. 
Just like in the simulations of the pulsed system, the Floquet operator is calculated on a truncated basis and some 
caution has to be exerted to avoid spurious results. A study of the quasienergy curves themselves can help us in 
deciding what parts of what curves to trust [24]: see Appendix A. 

2. Quantum nonlinear resonances 

The behaviour of a grouping of levels strongly interacting (repelling each other) at or close to zero microwave field 
is particularly evident in Fig. 1. I have marked them as darker lines. If we look at the (zero microwave field) energies 
of these levels we see that each level is close to a quantum resonance 

AE = E n+r - E n = sivh (4) 

with any of the other levels of the same grouping; the ratio s/r being equal to 1 for any two levels of the grouping. 
This grouping represents on the (E, F) plane what Ref. [25] calls a nonlinear quantum resonance: a finite (because 
of level anharmonicity) number of quantum states whose interaction is due to the existence of a (primary) classical 
resonance zone, in this case the u>' = 1 resonance zone. The structure these levels form can be described as an 
excitation of the system induced by the perturbation: Ref. [26] calls it a "Floton" . The state of the Floton is decribed 
by the new quantum number k given in Fig. 1 next to the free atom one n; a simple interpretation of it in terms 
of pendulum approximation is given in Appendix B2. The most stable state is the ground one (k = 0); the excited 
states k = 1,2, .. are progressively shorter lived [27], as long as we consider only the properly resonant states: those 
for which the inflection point of the corresponding quasienergy curve in Fig. 1 lies to the left of the considered value 
of the microwave field strength F (see Appendix B2). 

3. Husimi function 

One of the most widely used representations of a wavefunction in phase space is the Husimi function [28] . Given 
a quantum function its Husimi function at a point {(p) , (q)} in phase space is the projection (in atomic units) 
Ph((p) , (<?) , t) = | ( ( P{p){q)\' l P('t)) 1 2 °f orL a (minimum uncertainty) packet (j>{p)(q) centered on that same point 
in phase space and called the "coarse graining function" . In action-angle space, due to the periodicity in the angle 
variable, the standard (position-momentum space) choice of coarse graining function [28] is not possible. Following 
Ref. [29] I therefore take 

(an+l) 

^ = S "= o[ 2.r(an + l) (/> 6 ] 6 ^ 
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Tests with different values of the width parameter a have given - both for Ref. [29] and for me - Husimi functions 
having essentially the same shape when a < n . For n 3> 1 the classical zones in phase space are reflected in the 
character of the Floquet eigenstates (FE), especially as seen through their Husimi functions (see refs. [24, 29] and - for 
semiclassical FE - Ref. [8]). Thus regular FE can be either resonance island FE or modulated atom FE. Irregular FE 
include the separatrix FE, possibly "scarred" by a high electron probability density in phase space along an unstable 
periodic orbit [12, 27]; in sectioned phase space this corresponds to a maximum of the Husimi function on the unstable 
fixed point which (surprisingly) has a simple interpretation in pendulum approximation [26] . At intermediate values 
of microwave field strength, these scarred separatrix FE exhibit partial quantum localization along the unstable 
periodic orbits (rather than full localization on nested vortex tubes), and otherwise are distributed within a bounded 
region of phase space classically occupied by chaotic trajectories [8, 27]. The comparison of Husimi functions for the 
quasienergy states with the classical surfaces of section at the same values of u, F and <f) has often proven quite 
instructive. Comparison of Husimi functions and classical instantaneous ensembles for the pulsed simulations is also 
possible, but I found them less useful in the present study. 

III. RESULTS 

In Ref. [13], the large secondary resonance island chains encircling the stable fixed point of the Poincare map at 
the center of the lo' = 1 resonance zone were identified as the cause of the sequence of steps in the classical ionization 
probability vs. peak microwave field F™ ax graphs, of which three examples are shown as full lines in Fig. 2. These 
chains are born from the distruction of those periodic orbits of the Poincare map laying within the uj' = 1 primary 
resonance zone which have periods whose ratio to the microwave period equals the ratio q/p of two (relatively prime) 
integers [13]. As Fq is increased, The most visible chains have p = 1 ("integer resonances") and the number q of 
individual islands in them decreases with F [13]. The size of an island chain grows rapidly with F above a bifurcation 
point F 0q / p in qualitative agreement with the simple pendulum approximation model described in appendix B 1. 

The arrows in Fig. 2 show the field strength values at which the p = 1 island chains destabilize (become completely 
chaotic) and are destroyed, for q = 3 to 6. Each of these values is very close to a step onset, strongly suggesting a 
connection between the ionization steps and the secondary island chains; the connection is confermed by the evolution 
in phase space of the classical trajectory ensemble. Four examples - of the many I have calculated in Ref. [16]- are 
given in Figs. 3-6: in all cases we are at the base of an ionization step and a sizable part of the ensemble (marked by 
bigger dots) finds itself at the peak of the pulse inside the islands of some secondary chain: q = 6 for Fig. 3, q = 5 
for Fig. 4, q = 4 for Fig. 5 and q = 3 for Fig. 6. Together with the big whorl of the points trapped within the 
primary island, the marked points make up most of the sub-ensemble surviving the microwave pulse, as can be seen 
from the last snapshot of each sequence. Whorls and tendrils developed by the ensemble at intermediate times during 
the rise of the pulse indicate the influence of other secondary chains with higher g's. The relevant chains were in 
each case identified by comparing the snapshot of the ensemble with the constant field stroboscopic surface of section 
taken at the same value of Fq; examples are shown in Figs. 7-10. These comparisons reveal the following classical 
time sequence of events responsible for a step in Figure 2 (a more schematic description was given in Ref. [13]). The 
ensemble initially is located on a nonresonant phase-space KAM curve which spans the full range [0, 27r] of the angle 
variable. The pulse peak field F™ ax must of course be high enough that by the time we are at the peak of the pulse 
the locally chaotic region of the principal primary resonance has merged with the global chaos above (the last KAM 
curve separating them has broken) so that classical ionization is possible. 

A) For F (t) near 0.01, the ensemble almost entirely crosses the locally chaotic region to enter the uj' q = 1 island 
region; the wide "nose" the ensemble develops at the crossing (see the second snapshots in Figures 3-6) tells us that 
after the crossing the spread in resonance action J (action of the slow motion around the stable fixed point of the 
Poincare map) is quite wide. 

B) After the crossing there is a time interval of slow near adiabatic ensemble evolution, having an approximate 
frequency given by eq. (B6). The ensemble is at the same time partially "adiabatically squeezed" in the angle variable 
[30] toward the center of the island (fourth to sixth snapshots in Figures 3 and 4, third to fifth snapshots in Figure 
5 and third to fourth snapshots in Figure 6). If \ui' — 1| < 0.05, much of the ensemble approaches and stays close to 
the center of the island. 

C) A bifurcation sequence with decreasing q has been occurring with increasing Fo(i). For each q the size of the 
corresponding island chain grows rapidly [13]. Within a field strength increment between 5% and 20%, this size 
reaches a value about equal to the nominal size of the entire primary island region (An = 8 to 11 depending on F (t)). 
The size of the individual islands of the chain grows too; for a given size AI of the chain the size of its islands appears 
to increase with decreasing q. 

D) A study of surfaces of section at F < 0.01 (crossing field strength) shows that the secondary island chains with 
q > 8 are all born before the crossing; most of them (those with q > 12) have moreover already grown to full size and 



6 



destabilized. The remaining secondary island chains can and - for q < 8 - do influence the evolution of the ensemble 
as can be seen in the superpositions of ensembles over instantaneous surfaces of section shown in Figs. 7-10: parts 
of the ensemble are "trapped" by these chains. By "trapped" I mean that once the points of a part of the ensemble 
are inside the islands of a chain, they remain inside those islands and follow the chain in its motion (during the rise 
of the pulse) toward the periphery of the primary resonance zone [31]. Note that the direction of this motion is the 
opposite of that of the adiabatic "squeezing" we observed in step B. This inversion of motion of parts of the ensemble 
within the primary resonance island is a clear indication of their trapping by the secondary resonance zone. Finally, 
since the span, in the resonance action J, of the islands of the secondary chain is usually not big, it is important that 
a sizable portion of the ensemble be at some time within the same small range of J so as to be able to be trapped by 
the islands of the secondary chain. This is made possible by the "squeezing" we have seen at point B. 

For each of the four cases in Figs. 3-6 I have marked with darker dots that part of the ensemble which at the peak 
of the pulse is within the relevant island chain: the motion of the "trapped" part of the ensemble toward the periphery 
of the primary resonance zone is particularly clear in Figures 5 and 6. It is also evident that different island chains act 
on different parts of the ensemble: the higher is the q of a chain the further from the center of the primary resonance 
zone is the part of the ensemble it "traps" . Parts of the ensemble closer to the center of the resonance zone do not 
"see" the chain when it grows through them: they are traversed diabatically (motion within the secondary islands is 
too slow to be "seen" by the points of the ensemble for the chosen switch-on time). Only at a certain distance from 
the center of the primary resonance zone a given island chain becomes "visible" to the ensemble, and this leads us to 
believe that, together with the size of the chain and the size of the islands of the chain, also the (average) frequency 
of motion within those islands must increase with Fq (just like the frequency of motion within the primary island: 
see eq. (B5)). Moreover, for a given pulse length, the distance from the center at which this happens increases with 
the q of the chain, suggesting that, for any given distance from the center, the higher is q the lower is the frequency 
of motion on the islands of the chain. 

Now that we have part of the ensemble trapped in the secondary chain, we have two possibilities : E and E'+F'. 

E) As long as F™ ax remains smaller than the one necessary to destroy the islands of the secondary chain, the only 
part of the ensemble which can and does ionize is the one in the chaotic region of the primary resonance zone; no 
further ionization is possible. We are on the plateau of the step. 

E') If we instead further increase F™ ax , the island chain becomes locally chaotic, and then "globally" chaotic (that 
is: it merges with the chaotic band surrounding the primary resonance island and finally with the global chaos above), 
all within a further field strength increment of about 5%: rather abruptly a sizable part of the phase space becomes 
globally chaotic. It is this sudden increase of the globally chaotic region that causes the step in ionization probability. 
The growth of the band of local chaos is reflected in the changed appearance of the parts of the ensemble trapped in 
the chain: whorls, connected to motion within the stability islands, are substituted by tendrils, connected to motion 
in the chaotic band, the last two snapshots in Fig. 7 show such a transition for the q = 6 island chain; the association 
of tendrils with the destruction of chains is also clear in the second snapshot in Fig. 9 (again the chain is the q = 6 
one) and in the third (q = 5), fourth (q = 4) and sixth (q = 3) snapshot in Fig. 10. 

F') The chaotic trajectory subensemble released by the island chain destruction event E' produces ionization during 
a time interval of about 30% of the pulse length, within the relatively constant central part of the half sinewave pulse. 

The examples I have given in Figures 3-6 and 7-10 make clear that "trapping" in the secondary island chains is 
possible for the experimental pulse times of about 131 microwave periods. Longer pulse times mean slower growth 
of F (t) resulting in even more effective trapping, as can be seen from Fig. 14 in Ref. [32]: the sharpness of the 
ionization probability steps improves with increasing pulse times. 

An important implication of the above scenario is that, for fixed values of the static field F* (rescaled at the resonant 
action I r defined by a; 7^/(1 — 3F S I^) = 1), the position in F™ ax of the steps depends only on the destabilization of 
the island chain and not on the initial conditions of the ensemble. The initial action (and the pulse time) on the other 
hand determine the heights of the steps. As an example I have plotted classical ionization probability data for rescaled 
frequencies ui' between 0.9325 and 1.0666 first as a function of F™ ax (Fig. 11. a) and then as a function of F™ ax (Fig. 
11. b): notwithstanding the great difference in height, the position of the steps in Fig. 11. b is approximately the same 
(F^ is slightly different from curve to curve); while this is clearly not the case in Fig. 11. a [33]. 

Parallel quantum numerical simulations are shown as dashed curves in Fig. 2. In Figure 2. a. The classical steps 
labeled q = 3 — 6 are reproduced quite well. They appear even sharper in the quantum simulations than they are 
in the classical ones. On the other hand the oscillations of the plateau just before the q = 3 step are much more 
pronounced than the classical ones. The absence of a q — 2 quantum step in the figure most likely reflects the fact 
that the q = 2 island chain in (sectioned) phase space is not large enough to significantly support a Floquet eigenstate 
(even if their height is comparable to that of the islands of other chains, the two islands are extremely narrow: never 
more than about 0.1 radians in the angle variable). Figure 2.b shows another case where the above comments again 
apply: the quantum curve follows quite closely the classical one up to the q = 3 step but then directly rises to 100% 
ionization probability, skipping the q = 2 step. Figure 2.c finally shows the quantum and classical steps for a third 
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case; this time \ut' Q — 1| = 0.073 is big enough that the ensemble does not penetrate all the way to the island's center. 
I believe this to be the cause of the absence of the q = 2 step in the classical curve: by the time the q = 2 island 
chain is born, the primary island is so small that all of the ensemble is outside of it in the chaotic region and cannot 
be trapped. 

The good agreement between classical and quantum simulations is not limited to the final ionization probability; 
the evolutions during the pulse are remarkably similar too. In Fig. 12 I show a comparison of the quantum and 
classical ionization probability versus time: the shape of the curves is always very similar even if the final ionization 
probability is often not exactly the same. 



To explain this agreement between classical and quantum ionization results, we now explore the connection be- 
tween higher order classical resonances and quantum resonances at nonzero microwave field (avoided crossings in the 
quasienergy plots). We start with some general considerations which will then allow us to locate the avoided crossings 
related to the secondary resonance zones within the principal primary resonance island [16]. 

For simplicity, let us start considering an isolated avoided crossing of two energy levels E(n i} F), i = 1, 2 with F 
some perturbation parameter. We assume that we have been able to separate the complete Hamiltonian H of the 
system in two parts, Ho and Hi, such that only the interaction responsible for the separation of the levels at the 
crossing is included in Hi and the levels E {m,F), eigenvalues of H , actually cross at the field F c of closest approach 
of the levels E(rii,F) of H. We moreover assume that E (m,F) is the quantized version of a classical function 
Eo(I, F) which is smooth for / G [m, ri2\- Since the classical phase space can be divided into regions corresponding to 
different kinds of motion and discontinuities in action occur at the borders of these regions [34] , the above condition 
requires us to assume that for all relevant F's the two levels of H we are considering correspond to the same kind 
of classical motion (their support is in the same region of the classical phase space of Ho). When these conditions 
are verified, at the field strength value F c where the levels Eo(rii, F) cross the average slope of the function E (I, F c ) 
for / G [jii,n2] is 0. There must therefore exist an action I r G [^1,712] verifying the classical resonance condition 
dEo(I,F c )/dI\i r = (Rolle's theorem). For such an avoided crossing and for small values of the perturbation, a 
relationship has been derived in Ref. [35] expressing the splitting at the avoided crossing as a function of the area of 
the related resonance zone, both functions of the interaction Hamiltonian Hi. 

If the perturbation Hamiltonian Hi is time dependent and periodic with frequency u>, then the avoided crossing is 
between quasienergy states, the quantum resonance condition for the (no more crossing) energy states of H is eq. 
(4) and the classical resonance condition is dE (I , F c )/dl\i r — suj/r. 

Let us now take H to be the "pendulum approximation" Hamiltonian 



which describes the resonant motion within the principal primary resonance zone s/r = 1 (see appendix B2). Its 
numerically calculated levels for the parameter values of Fig. 2. a arc shown in Fig. 13 (again see appendix B2 for 
details). After each level's inflection point (which is when the level gets "trapped" [25] into the region of libration 
motion), the functional dependence of the energy of the eigenstates of H^ on the pendulum quantum number k is 
that of the classical energy on the action J within the (primary) "resonance" zone of libration motion (see fig 23 in 
appendix B2). Since this dependence is smooth we can apply the argument above: crossings between states with 
Afc = q and energy distance AE — pu; are related to the f2 = u>p/q classical secondary resonance condition where fl is 
the "slow motion frequency" \dH/dJ\. In Fig. 13 I have marked the crossings of the pendulum approximation levels 
with different symbols for the different p/q ratios (connected by lines to guide the eye); arrows mark the classical 
bifurcation fields eq. (B7) for the same p/q ratios. As expected from the argument above, all the crossings with a 
given p/q ratio lay to the right of the (pendulum approximation) bifurcation field for the secondary resonance zone 
with the same p/q ratio. Also, the crossings with the same p/q ratio are not all at the same field but the higher the 
average k of the levels crossing the higher the field at which the crossing takes place, in agreement with the moving 
(with increasing F) of the center of the classical secondary resonance zone toward higher pendulum actions [16]. 

We can now turn to the quasienergy levels of the full problem, of which the pendulum ones are a good approximation 
near the s/r = 1 quantum nonlinear resonance at low microwave field strengths [36]. The true levels' avoided crossings 
corresponding to the approximate levels' crossings, even if shifted (at high microwave field strengths) with respect to 
these latter ones, are related to the same secondary resonance zones and therefore lay to the right of the numerical 
classical bifurcation fields for the same p/q ratios, as can be seen in Fig. 1 (where we use the same conventions as in 
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Fig. 13) [16]. 
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The identification of the crossings related to these secondary island chains is supported by the Husimi functions 
I calculated for the Floquet eigenstates undergoing the avoided crossings: sonic of them are shown in Figs. 14- 
18. Before discussing them let us first introduce one last concept: that of crossing region. For F -C F c the two 
quantum states E(n,i,F) are practically noninteracting and coincide with E (n i} F)-. approaching F c they begin to 
interact more. This interaction is reflected in the mixing of the two states which is maximum at F = F c where the 
two eigenstates of H undergoing the avoided crossing are linear combinations in equal parts of the two (crossing) 
eigenstates of H n . After that, the mixing decreases until, for F 3> F c , they are again the two pure E (n i} F) states 
we started with but now the state m has the spatial structure n 2 had before the crossing and viceversa, as if the two 
levels had switched. The quantum resonance itself as defined by cq. (4) is therefore only the central point of a region 
in F where the two levels interact with each other. The values of F such that the separation of the noninteracting 
levels E (ni, F) equals the interaction term can be taken as approximate boundaries in F of the region where there 
is significant mixing [24]; we call this region the crossing region. Away from crossing regions the Husimi functions 
of the primary resonance eigenstates appear approximately as circular ridges, apart from the k = one which is 
approximately a Gaussian peak [24] . All these ridges have their center in the center of the classical principal primary 
island and the higher the resonance quantum number fc, the higher their average radius. Examples are the Husimi 
functions before and after the avoided crossings shown in Figures 14-18. At a crossing related to classical secondary 
chain of q islands, the exchange of the spatial structures of the two eigenstates involved happens in a characteristic 
way: the outer ridge moves inwards avoiding ("seeping" around) the q regions of the islands of the chain while the 
inner ridge moves outwards right through these regions, avoiding the regions around the unstable periodic orbit of 
the island chain and sometimes developing q peaks corresponding to the islands of the secondary chain. Often this 
q-fold structure of the states is not equally visible for both states at the same value of the microwave field strength 
F: for most of the cases studied the structure is most visible for the outer states at slightly lower values of F than 
for the inner state. Fig. 14 is a particularly clear case: the Husimi functions of the fci = 2 and k 2 — 7 resonance 
states both develop at their crossing the expected fivefold structure. Figure 15 shows instead the (q — 7) crossing of 
the ki = 3 and k 2 — 9 states: since the k 2 = 9 state is here almost a separatrix state (its Husimi function has not yet 
taken the typical resonance state shape and is peaked around the unstable fixed point of the Poincare map [27]) we 
are nominally outside the range of applicability of my semiclassical argument, still the sevenfold structure is visible. 
In other cases the process can be complicated by the interaction with other states but the g-fold structure is always 
present. Figure 16 shows the (q — 6) crossing of the states fci = 2 and fc 2 = 8 where the lower "adiabatic" state (with 
respect to the q = 6 crossing) is made up of the two interacting states B and C. In Fig. 17 (q = 5, k\ = 0, k 2 = 5) 
it is the upper "adiabatic" state which instead is made up of the two states A and C. Finally Fig. 18 shows a case 
(q = 4, ki = 0, k 2 = 4) where the three level interaction results in the adiabatic state "B" assuming before the crossing 
the role of upper "adiabatic" state and after the crossing the role of lower "adiabatic" state, thus maintaining the 
same resonance quantization k 2 = 4. These Husimi functions make clear that at the crossing regions I have marked in 
Fig. 1 the crossing eigenfunctions do see (and reflect in their structure) the q dips of the classical dynamical potential 
[37] evidentiated in the surfaces of section by the q islands of the secondary chain they are related to. 

We now have to show that these crossings are in effect those where most of the population transfer causing the 
ionization steps happens during the microwave pulse. To do this we plot the projections of the wavefunction on 
the instantaneous Floquet eigenstates at each period of the microwave during the pulse. For the four values of the 
peak microwave field strengths of Figs. 3-6 these projections are shown in Figs. 19-22. They give us a dramatic 
representation of the importance of the avoided crossings we have individuated (marked here with the same symbols 
as in Fig. 1). After the initial spread of the population over the states k = — 3, the population is further spread 
out to high k short lived states through these avoided crossings or sometimes through avoided crossings which can 
be related to secondary island chains with higher q than those marked. So in Fig. 19 the two (couples of) crossings 
which are not marked are those between k — 1 and k — 9 and between k — 2 and k — 10; both can be reconducted 
to the q = 8 secondary island chain but only by analogy with the cases studied, since the crossings happen early 
in the k = 9 and k — 10 curves and their WKB quantization is not yet the resonance one (i.e. depending on the 
quantum number fc). The (q = 8) crossing between k = 2 and fc = 10 appears again as significant in Fig. 20 but in 
Fig. 21 the transfer of probability at that crossing has dwindled to almost nothing and the only significant crossings 
are the marked ones. Finally in Fig. 22, corresponding to a case where the classical principal primary resonance zone 
is almost completely chaotic at the peak of the pulse, the marked crossings are again the first relevant crossings and 
the only ones to cause significant transfer of probability; but several other crossings cause some (minor) transfer of 
probability at the peak of the pulse. The most important of these latter crossings are related to classical secondary 
"fractionary" resonances: for example the crossings marked by crosses (from top to bottom, between fc = and fc = 9, 
between k = 1 and fc = 10 and between k = 1 and fc = 11) have p = 2. There also appears to be some overlap of 
significant crossings each involving the same level but related to different classical resonance zones (for example the 
crossings between fc = 1 and fc = 5 and between fc = 1 and fc = 11 in the bottom right corner of the figure); this 
is by some authors thought to be a sign of an analogous overlap of the related classical resonances and therefore an 
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expression in quantum dynamics of local classical chaos [38] . To emphasize the step-like character of the evolution of 
the k = — 3 states, I give in Table I their population loss at each of the relevant crossings. 

IV. DISCUSSION 

We have thus shown the relationship between some strictly quantum processes (avoided crossings) and classical 
structures; but why do isolated avoided crossings so well simulate the classical behaviour in the cases we studied? To 
understand this we should try and get an idea of the circumstances under which the classical and quantum dynamics 
can be expected to result in similar evolutions of the observables. 

Let us start considering the system at a given microwave field strength; when can we say that the properties of the 
quantum system reflect the classical resonance structure? Away from any crossing region related to a given secondary 
classical resonance, the quantum system is by definition not affected by that classical resonance zone. Only for values 
of F falling within one such crossing region the quantum eigenstates reflect the classical resonance region, as we have 
seen. If the crossing regions related to the classical resonance zone are well spaced, it will therefore be only by a 
chance choice of F that the quantum dynamics will be able to "feel" the classical resonance zone. Only when the 
crossing regions related to the classical resonance zone overlap this element of chance is lost. 

If we look at Figure 1 we see that there is hardly any overlap of the crossing regions for any of the p/q ratios marked 
in the figure; unless F is carefully chosen we therefore expect the fixed microwave field strength quantum dynamics 
not to be affected by the classical secondary resonance zones for our choice of no; in other terms: we cannot be sure 
of having any quantum state reflecting (mimicking) the classical secondary resonance structures. This is essentially 
[39] the base of the objections which have been raised to a classical interpretation of the steps [32], notwithstanding 
the good agreement of classical and quantum ionization probabilities. I disagree on the ground that arguments for 
an atom interacting with a constant amplitude microwave field are of limited relevance to the present pulsed system: 
since the microwave field strength F(t) changes during the pulse, the system in general meets not just one but many 
avoided crossings related to a single classical resonance (see the examples in Figs. 19-22). Most of these crossings 
are between completely different states (e.g., for the q = 6 classical resonance, between the states k = and k = 6 , 
between the k = 1 and the k = 7 ones and so on ...) and each of these avoided crossings is a potentially very quantum 
event as we do know that, when a crossing is met by the system twice (once during the rise and the second time during 
the fall of the pulse), the different phases accumulated by the two parts of the wavefunction (which split at the first 
crossing) cause interference at the second crossing, resulting in Stueckelberg-like oscillations in the population of the 
quasienergy states as a function of F max [40]. We do not see any trace of this quantum behaviour. I believe the cause 
of this to be twofold. First of all, for the current value of the pulse length T, the separatrix crossing event early in the 
pulse almost immediately significantly spreads out the population on states with resonance quantum number up to 
k = 3 as can be seen in Figures 19-22 [41]. This means that, even if there were Stueckelberg oscillations, the ionization 
signal vs. F max would be the sum of oscillations due to the couples of crossings met by all four these states, each of 
the four states meeting a sequence of them. The period of each of these oscillations depends on the phase difference 
accumulated by the two crossing states between the crossings, the amplitude on the crossing width (and the speed of 
the crossing) and the phase on the F max when the crossing is first met by the system. That these oscillations should 
sum up to a coherent oscillation, rather than average out would be a rare event indeed. Were this the case in our 
simulations, we should be able to see some sign of these oscillations in the populations of the single states. Namely 
we would expect to sometimes see at the avoided crossings in the falling edge of the pulse a return of probability from 
the high k states to the four initially populated ones (k = 0, 1,2 and 3). While this transfer is visible at the very 
end of the pulse among the four states themselves, we have no evidence of this happening at the avoided crossings we 
are considering, as can be seen from Table I. A posible cause (our second reason for quasi-classical behaviour) is the 
very low population on the high k states at crossings on the falling edge of the pulse: at low peak microwave field 
strengths only high q (narrow) avoided crossings are met so that the population transfer on the rising edge of the 
pulse is very small. With increasing peak microwave field strength wider crossings, corresponding to lower g's, are met 
and population transfer on the rising edge becomes substantial, but so does also the decay rate of the high k states 
between the two crossings. The importance of the spread of the initial population can on the other hand be inferred 
from the comparison made in Ref. [32] between the classical and quantum final ionization probabilities for the two 
pulse times T — 131 (our case) and T = 300 microwave periods. In both cases the agreement increases with increasing 
principal quantum number n ; but a good agreement is reached at lower values of no for T = 131. This latter is the 
case where we expect more primary resonance quantum states to be excited at the crossing of the principal primary 
island separatrix. In this picture, rather than being traces of Stueckelberg oscillations, the fluctuations of the plateaus 
of the quantum steps in Fig. 2 are likely due to the variations of population transfer at the avoided crossings during 
the rising edge of the pulse. For any given q, the first crossing to be met involves the k = state and, when first met, 
happens at the peak of the pulse, so that the interaction time is long and the population transfer high [40]: we have 
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a ionization peak. With increasing peak microwave field strength the interaction time at the crossing decreases and 
so does the population transfer; but then the k = 1 state has a crossing with the same q and a new ionization peak 
appears on the tail of the k = one, and so on. 

V. CONCLUSIONS 

My quantum simulations show that each of the experimentally observed ionization steps is due to a group of avoided 
crossings between couples of Floquet states belonging to the same quantum nonlinear resonance. These crossings are 
in their turn related to the classical secondary resonance island chains surrounding the principal primary resonance, as 
the quantum and classical resonance conditions are formally the same. Finally, the above island chains are responsible 
for the steps observed in the classical simulations, as their distruction releases into the globally chaotic phase space 
region (and therefore to eventual ionization) those parts of the electron ensemble previously trapped within them. 

The good agreement between quantum and classical simulations is a consequence of the spread - early during the 
pulse - of the quantum probability to several resonant Floquet states. This spread allows more than one of the avoided 
crossings related to each of the secondary island chains to redistribute the electron wavefunction to short lived states, 
thus contributing to the ionization step corresponding to that island chain. 

My study allows us to conclude that classical phase space structures too small to influence (in most cases) the 
quantum evolution for a constant amplitude microwave field, can -in a pulsed regime- generate interesting quantum 
dynamics. Constant amplitude arguments to determine their relevance should therefore be used with extreme care, 
as they can lead to wrong conclusions. 

VI. ACKNOWLEDGEMENTS 

The author wishes to thank C. Rovelli, G. Mantica and S. Locklin for useful comments and suggestions and the 
latter also for the use of his computer for the quantum numerical simulations presented in this paper. Special thanks 
to J.E. Bayfield as advisor of my Ph. D. thesis, amply quoted in the present paper. 

APPENDIX A 

I discuss here the choice of integration parameters for my quantum simulations [16]. Let us start with considerations 
about the extremes of the truncated basis I used. An occupation probability (at every time during the simulation) of 
the order of 10~ 10 on the lowest two or three states can be considered a sign that the basis lowest level n min is small 
enough. For all of my runs with n in the range 60 to 72, n min = 30 satisfied this condition. 

For the upper bound the matter is more delicate. The spectrum of atom in the static field alone consists of 
resonances (bands) whose number is not well determined, as increasingly wide ones exist with energies extending 
above the classical static field barrier. Even if the free atom states tp n of which my basis is composed extend to 
infinity their probability distribution is centered around (z) n — 3n 2 /2 with a spread Az n = n 2 /2. This means that 
truncating the basis puts us in a "fuzzy box" . The potential is therefore a kind of "double well" where not only the 
width but also the depth of the second well varies with the number of levels considered and the calculated spectrum 
shows couples of quasi degenerate double well states instead of bands. The narrow separation of the couples of levels 

1/2 

- up to just below the classical ionization thcshold E rnax = —2F S - tells us that in the static field alone there is very 
little coupling between the two sides of the well (as expected for an asymmetric double well [37, 42]). In my pulsed 
simulations, I take into account the loss of probability induced by the static field by the substitution in equation (3) 
of E n with (E n — iT n /2), where the decay factors r„ are given by the (3D) formula from Ref. [43]; this is a rather 
crude approximation: the decay factors are for the Stark states of the problem and I apply them to free atom states 
but still using Stark perturbative energies in calculating them. This would place us in trouble with states with Stark 
energies around the maximum of the atom in static field potential E max , were it not for the sharp increase there of the 
decay rates from practically zero to practically infinity over a range of two or three states. This makes me confident 
that the errors connected with the use of the wrong states will concern only a few states close to the ionization border 
which never in our simulations get significantly populated. The continuum above the static field ionization threshold 
gives no cause for worry: for most of my simulations the uppermost level is n max — 221. This gives us about 60 states 
above the static field ionization threshold, extending in energy about 5.6 * 10~ 5 a.u. above that threshold. For the 
frequencies of my simulations (u) < 2.7 * 10~ 6 a.u.) this span in energy is certainly sufficient; moreover, due to the 
big decay rates of the states of this "pseudo-continuum" , the system will not be able to resolve the individual states. 
More serious cause of worry is that my procedure actually cuts away pieces of the "double well" states and this might 
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cause interference effects altering the final result of the simulation. On the other hand such interference effects should 
be highly sensitive to changes of the basis and extensive tests show this not to be the case. 

When the microwave field is switched on, some indication on the validity of the above approximation can be 
extracted from the quasienergy curves. Ref. [24] has shown that the quasienergies calculated on a truncated basis are 
subject to the unphysical constraint that their sum is a constant independent from F but dependent on the number 
of states in the basis. Whatever our choice of basis, certain curves (at the top of the basis) are going to be completely 
wrong. Now, with increasing microwave field some of the quasienergy states move up and some down; the wrong 
states are therefore bound to undergo a series of avoided crossings. Since the coupling among the states increases 
with the microwave field strength, the wrong curves soon or later interact (have a pronounced avoided crossing) with 
other curves. The avoided crossings of these other curves for fields above the field at which this interaction takes 
place therefore become in their turn unreliable and so on, to the point that all the quasienergies of interest bear no 
resemblance to those of the complete problem [44]. Often it is still possible to extract some information from the 
"wrong" sections of the curves, but it will be mostly qualitative. There are some rules of the thumb to recognize the 
completely wrong curves. Levels close to the top of the basis are always suspect. Moreover when the nearest neighbour 
free atom (or, in our case, atom in static field) energy spacing is more than ivh one expects [24] the quasienergy curves 
to bend down and viceversa. Curves which do not follow this pattern are usually suspect. The most reliable way of 
finding the wrong curves is still, on the other hand, to perform test calculations on different bases. For my calculations 
the free atom basis used was the same I used for the pulsed simulations. For frequencies of interest the lowest level 
n min = 30 (and those just above it) has no wide avoided crossings and is practically at the same quasienergy over 
the entire range of field strengths spanned by our simulations. Tests where I included all the states down to n=l 
show that all levels with lower quantum numbers behave the same way. Tests where I reduced the upmost level from 
n ma x = 221 to n max = 98 show that, for the microwave frequencies of interest here, the behaviour of the levels near n 
are very similar over the range of field strengths spanned by my simulations. On the other hand the avoided crossings 
of the relevant states with the wrong curves at the top of the basis are very narrow for the bigger basis but become 
much wider for the smaller basis. 

We have seen that the free atom basis somehow functions as a Sturmian one does in absence of a static field, 
reproducing the resonances above the classical static field ionization threshold; but we have also seen that truncation 
of this basis means the presence of an artificial second potential well with unphysical finite number of discrete states 
in it. One should therefore wonder about the impact of these levels on the levels of interest. Curves corresponding 
to states deep in the second "well" tend to have very little interaction with any other state and higher second "well" 
states tend to have significant avoided crossings only with each other (avoided crossings wide enough that for our 
pulse lengths they are not crossed diabatically and some population exchange takes place [40]). For the range of 
microwave field strengths of interest, the only "second well" states having (significant) avoided crossings influencing 
the "bound states" we are interested in are those close to the classical static field ionization threshold. 

Another parameter which can create serious problems is the integration step At: it must be small enough so that 
the spurious frequencies 2irk/At induced by the time discretization [15] do not cause transitions between the basis 
states. Approximating the basis upper limit with the free atom continuum, this means that any At < 47rn^ lin is 
sufficient to inhibit spurious transitions even between the very little populated extremes of the basis. For n min = 30, 
the above condition is At < 10 4 a.u.. For most of my simulations At = 5 • 10 3 a.u.. 

APPENDIX B 

1. An evaluation of the birth and growth of the secondary island chains using the classical pendulum 
approximation for the principal primary resonance island [16] 

When F(t) is a constant, one procedure often used to study the Hamiltonian (1) is to locally approximate it with 
an integrable one around any resonant action I r defined by w 7^/(1 — ZF s lf) = s where s is an integer [17]. The 
resulting Hamiltonian is that of a pendulum describing the slow motion inside the resonance island; for the principal 
primary resonance s = 1, it reads 

=-(^+ a cosv (Bl) 

where 

/3-3(l + F 5 / r 4 ), a = 0.325F , 
ip = 9 — ujt, 
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and 

The action angle variables of the Hamiltonian (Bl) itself are well known. For rotating motion outside the separatrix 
the (rescaled) action and angle are [45] 

j_ 4 fa\ 1/2 E(i?) 



ir\0j R ' 



where 

R= ' 



1-H/aJ |sin($/2)|' 

$ being the (half) amplitude in <p of the oscillations. E and K are complete elliptic integrals [45], F is the incomplete 
elliptic integral corresponding to K, and the + sign refers to i > 0. Analogously - for librating motion inside the 
separatrix - [45] 

8 fa\ 1/2 r (R 2 -l)K{l/R) 



E(l/R) 



R 2 



RF(cp/2,R) 

^- J = ±7r 2K(l/fl) (B3) 
The frequency of motion in this latter case, known as "slow motion frequency" [46], is: 

= = \dH^/dJ\ = 7T (a/3) 1/2 /2K(1/ J R) (B4) 

which is maximum for l/R = where it becomes the usual "small oscillations" frequency 

n max = (a/?) 1/2 . (B5) 

Since there is for every F a maximum slow frequency Vi max , the resonance condition for a secondary resonance with 
rotation number p/q can be verified only if 

n max > pcu/q = p(l - 3F 3 I?)/qll (B6) 

The above condition for the existence of a p/q secondary resonance zone can be written in the form of a condition on 
the microwave field strength: 

r> vlq r ~ 0.975g 2 (l + F s / r 4 ) ~ 0.975g 2 1 ' 

For F < F p / q the secondary resonance zone does not exist and it appears for F = F p / q at the center of the primary 
resonance zone, £l ma x being the pendulum libration frequency there. Contrary to the primary resonance zones, it is 
clear that since f2 = Q(J, F), the secondary resonance zones do move in phase space. To see how they move we can 
combine (for F > F p / q ) the p/q resonance condition (B6) with eq. (B4) expanded in powers of l/R = |sm($/2)|; this 
gives us the distance in (p of the secondary island chain from the primary fixed point: 

$ - (8AF/F p/q y/ 2 (B8) 

where AF = F — F p / q . We thus see that for every p/q the resonance starts for F = F p / q - called the p/q bifurcation 
field - at the fixed point itself ($ = 0) and then moves outwards toward higher $'s with increasing F. 

Eqs. (B7) and (B8) have been obtained in pendulum approximation, we therefore expect them to work only up to 
moderate fields; their failure for low q integer resonances - which appear only at high microwave fields - is particularly 
evident for eq. (B8) where the exponent 1/2 does not fit the numerical results for q < 4 [13]. Another evaluation 
of the bifurcation fields F p / q can be obtained [16] from a stability analysis of the Kepler map: an approximate map 
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for the classical dynamics derived in Ref. [17]. In the map the fixed point at the center of the principal primary 
resonance becomes unstable after the p/q = 1/2 bifurcation point. This is confirmed by my numerical simulations 
for the complete system which show no trace of the primary stable fixed point above the -F\/2 bifurcation point. The 
fractional deviation of the Kepler map bifurcation points from the numerical ones is approximately constant: they 
are always 25% below the numerical ones, while the pendulum approximation bifurcation points are always above the 
numerical ones and their deviation increases from about 5% for q = 7 and 6 to about 128% for q = 2. The introduction 
of a static field F$ = 8V/cm reduces in our case the bifurcation points by about 20%. Again, in all cases the results 
in pendulum approximation arc higher than the numerical ones and the results from the Kepler map lower; but while 
the deviation of the pendulum approximation results increases from about 15% at q = 6 to about 150% at q = 2, the 
deviation of Kepler map results from the numerical ones remains about 15% from q = 6 to q = 2. Since in absence of 
a static field the deviation was about 25%; this suggest that the linear approximation of the static field I used is an 
over-evaluation. 



2. quantum pendulum approximation and WKB quantization [16] 

Following Ref. [25] and [36], the quasienergies of the Hamiltonian (1) can be approximated by 



-(i) _ 



= (E R + W k )(mod toh). 



(B9) 



where Er is the eigenvalue of the (classically) resonant state \R > of the Hamiltonian for the atom in the static field 
alone; and Wk denote the eigenvalues of the matrix Hn) m defined as follows: 



H£l = [E n -E R -w{n-R)], 

L n,n+1 



tfl'L, =F(R\z\R+l) /2 



(BIO) 



While the diagonalization of the Floquet operator gives us no natural ordering of the quasienergies, the energies Wk can 
be ordered. A descending order in energy with k starting from zero allows us to identify the index k as the quantized 
version of the pendulum action J. Moreover, since eq. (B9) gives energies folded into the first Brillouin zone, the level 
crossings appearing as a consequence of this folding are actual crossings. Eq. (B9) has been shown to give for quantum 
numbers close to R quasienergy curves remarkably similar to those of the full problem [36] for n = 60 and rescaled 
fields up to about Fo = 0.04. I have numerically calculated the "quasienergies" (B9) for the parameters of Fig. 1; the 

matrix elements Hnjn between the Stark states were calculated by projecting the Stark states on the free atom basis 
n € [33, 94]. The result is shown as the dotted curves in Fig. 23. For comparison, I also plotted - as big dots connected 
by dark lines - the result of a WKB quantization of easily calculated approximations of the constant action Jo curves 
of energy versus microwave field strength F implicitly defined by eqs. (B3) and (B2). The approximations made were 
the following: The zero action curve is the energy of the fixed point: H = a. Near J = we can expand K(l/i?) and 
E(l/i?) in powers oil/ R (to the fourth order [45]) and then solve for H the resulting approximation of eq. (B3). For 
a given action J this approximation will get better and better for increasing a: in the high field limit equidistance 
in action will therefore mean equidistance in energy and the system will behave essentially like a harmonic oscillator 
with frequency £1 = (a/3) 1 / 2 . A similar procedure is used in the opposite limit (a — > 0) starting from eq. (B2). 
At the separatrix (H = —a), there is a discontinuity in J: J = 4(a//3) 1 / 2 /-7r outside and J = %(a/ (3) 1 / 2 /tt inside; 
but, from both sides, approaching the separatrix, dH/da tends to —1 so that the H{a) curves join smoothly (C 1 at 
least). Since their slope at that point is that of the separatrix itself, we also know we have an inflection point there. 
Finally the curve H(a) has a minimum given by the condition 2E(l/i?) = K(l/i?) which solved with the help of the 
numerical tables [47] gives us K(l/R) = 2.3211 and J = 8(a/(3)^ 2 (-H/a)E(l/R)/TT = 1.1128(3a//3) 1 / 2 . With the 
above approximations and using the WKB n-quantization outside the classical resonance zone and the fc-quantization 
inside, we have: 

Low field behaviour: 
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Minimum: 



High field behaviour: 
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In Fig. 23 the resonant action is 7 r = 69 and the zero microwave field parameters £fe(0) and |Afc(0)| = — Afe(0) = 
Wfe(0) = e r ^ £fe(0) used in the construction are given in table II (in atomic units). The agreement is quite good both 
at high and low fields; the shift of the high field curves increases with k but for each of these curves it appears to 
slowly decrease with increasing field. The inflection points, which classically are right at the separatrix, are at the 
center of the transition region between the two quantizations: to their left the states are essentially modulated Stark 
states while to their right they take the character of resonant states. In our semiclassical treatment we can therefore 
only say that they fall somewhere between the n-quantized point and the ^-quantized one (circled points connected 
by segments in Fig. 23) 



a. Semiclassical scaling of the quasienergy curves [16] 



Reference [48] suggests an approximate semiclassical scaling law: 

e n (F)/u,*' 3 * f^iF/w*/*) (B12) 

where the function / depends only on classically rescaled quantities. Eq. (B12) connects quasienergy curves at 
different frequencies: to get the same function f nu i/3 for a different quantum number n, lo has to be changed. The 
agreement of my classical pendulum approximation evaluations with the quantum curves confirms for the cases studied 
this hypothesis, but only up to a certain point: while my evaluation for the low field behaviour, the top equation in 
(Bll), scales exactly according to eq. (B12), the evaluations for the minimum and the high field behaviour depend 
not only on n but also on (k + l/2)/n, as can be seen from the two bottom eqs. in (Bll). 
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TABLE I: The change of population on the four central resonance states in Figures 19-22 at the relevant avoided crossings. 


The populations are in percent of the total unit population. The q for each crossing is given in parenthesis, an asterisc follows 
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TABLE II: Quantum parameters for the semiclassical curves in Fig. 23. (From Ref. [16]) 
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FIG. 1: Quasienergy curves for n = 17.50GHz and F s = 8V/ cm calculated on the basis n £ [30, 221]. The horizontal axis is 
the microwave field strength rescaled for n = 69. The darker curves are those related to the s/r = 1 quantum resonance. The 
labels (n, k) give the free atom quantum number n and the resonance quantum number k. The groups of avoided crossings 
related to the most visible classical secondary island chains with rotation number 1/q are indicated by the following symbols: 
q — 3 circles; q — 4 down triangles; q = 5 squares; q — 6 diamonds; q — 7 up triangles. To aid the eye, lines are passed through 
the members of each group. Also, arrows mark the classical numerical bifurcation fields, labeled by q. (From Ref. [16]) 
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FIG. 2: Classical (dashed lines) and quantum (full lines) steps for a) u' = 0.9771, b) u>' = 0.9994 and c) u' = 0.927; 
Fq = 0.035266. Classical destabilization fields for secondary island chains are indicated by arrows; they are in all cases very 
close to the onset of a ionization step. (From Ref. [16]) 
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FIG. 3: Snapshots of the evolution of the classical ensemble for the parameter values lj' = 0.9771, Fq = 0.035266 and 
pmax _ q Q25 The ratio of Fo(t) to F™ ax is indicated in the top left corner of each snapshot. Bigger dots mark the part of 
the ensemble within the q — 6 island chain at the peak of the pulse. (From Ref. [16]) 
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FIG. 4: Snapshots of the evolution of the classical ensemble for the parameter values lj' = 0.9771, F(f = 0.035266 and 
pmax _ o The ratio of -Fo(i) to f™ 01 is indicated in the top left corner of each snapshot. Bigger dots mark the part of 
the ensemble within the q — 5 island chain at the peak of the pulse. (From Ref. [16]) 
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FIG. 5: Snapshots of the evolution of the classical ensemble for the parameter values lj' = 0.9771, F(f = 0.035266 and 
pmax _ q -phe ra ^j f Fo(t) to F™ ax is indicated in the top left corner of each snapshot. Bigger dots mark the part of 
the ensemble within the q — 4 island chain at the peak of the pulse. (From Ref. [16]) 
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FIG. 6: Snapshots of the evolution of the classical ensemble for the parameter values lj' = 0.9771, F(f = 0.035266 and 
pmax _ q_q61. The ratio of Fo(t) to F™ ax is indicated in the top left corner of each snapshot. Bigger dots mark the part of 
the ensemble within the q — 3 island chain at the peak of the pulse. (From Ref. [16]) 
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FIG. 7: Comparison of snapshots from Fig. 3 (F™ ax — 0.025) with the surfaces of section calculated at the values of the 
instantaneous i*b(t); Fo is indicated at the top of each graph, its ratio to Fg lax is given in parenthesis. In the first snapshot 
three islands of the q = 8 chain are deforming the tails of the ensemble; in the second one there are traces of tendrils of the 
almost completely destroyed q = 7 chain, while the q = 6 chain is beginning to be felt; in the last snapshot the q = 6 chain 
is almost completely destabilized, but still some tendrils and (on the left side) the remains of a whorl, are visible. (From Ref. 
[16]) 
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FIG. 8: Comparison of snapshots from Fig. 4 (F™ ax — 0.029) with the surfaces of section calculated at the values of the 
instantaneous Fo(t); Fo is indicated at the top of each graph, its ratio to F™ ax is given in parenthesis. In the second snapshot 
the q = 7 chain is beginning to be felt; in the third it is the q = 6 chain that is deforming the tails of the ensemble; and in the 
fourth the q = 5 chain. (From Ref. [16]) 



26 




FIG. 9: Comparison of snapshots from Fig. 5 [Fq 1 ^ — 0.041) with the surfaces of section calculated at the values of the 
instantaneous Fo(t); Fo is indicated at the top of each graph, its ratio to Fq 1 ^ is given in parenthesis. In the first snapshot 
both the q = 8 and q = 7 chains sligtly act on the enemble; but in the second the tendrils of the q = 6 chain are well developed; 
in the third the active chain is the q = 5 one; and in the fourth the q — 4 one. (From Ref. [16]) 
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FIG. 10: Comparison of snapshots from Fig. 6 (F™ ax — 0.061) with the surfaces of section calculated at the values of the 
instantaneous Fo(t); Fo is indicated at the top of each graph, its ratio to F™ ax is given in parenthesis. In the first and second 
snapshots the q = 7 and q = 6 chains respectively are the cause of the deformations of the tails of the primary island whorl; in 
the third snapshot the q — 4 chain begins to act but the tendrils of the q = 5 one are still visible; so in the fourth the tendrils 
of the q = 4 chain; the fifth and sixth snapshots see the creation and development of the q = 3 tendrils (no stability islands are 
present for q — 3). (From Ref. [16]) 
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FIG. 11: a) The classical probability for microwave ionization of no = 69 atoms plotted as a function of the rescaled peak 
microwave field strength F™ ax . There is no alignment of the ionization steps of the 25 curves shown (at values of uj'q equispatiated 
between lu' = 0.9325 and uj' — 1.0666). b) The same but plotted as a function of peak microwave field strength rescaled at the 
resonance F™ ax . With this choice of the horizontal scale the ionization steps are aligned. (From Ref. [16]) 



30 




FIG. 12: Comparison of quantum (full line) and classical (dashed line) time evolutions for various values of F™ ax (indicated 
next to each plot). F,f = 0.035266 and lo' = 0.9771. (From Ref. [16]) 
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FIG. 13: Pendulum approximation quasienergy curves for the case of Fig. 1. The groups of crossings related to the most 
visible classical secondary island chains are indicated by the same symbols as in Fig. 1. Here the arrows mark the classical 
approximate bifurcation fields eq (B7), labeled by q. (From Ref. [16]) 
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FIG. 14: Detail of Fig. 1 showing one of the avoided crossings marked in Fig. 1 and the Husimi functions of the (eigen-)states 
undergoing this same avoided crossing: it is a wide crossing where no other interacting level intervenes: the third q — 5 crossing 
(ki = 2, k<2 = 7). The (adiabatic) states are labeled by letters and numbers indicate the (diabatic) resonance quantization k. 
Bigger dots mark on the quasienergy curves the points corresponding to the Husimi functions displayed. Lighter shades of gray 
between the level curves mean higher values of the Husimi functions. The levels of the curves is the same for all plots; the 
highest peaks in some of them are, as a result, out of range and appear as white areas. 
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FIG. 15: Detail of Fig. 1 showing one of the avoided crossings marked in Fig. 1 and the Husimi functions of the (eigen-)states 
undergoing this same avoided crossing: it is a narrow crossing where no other level intervenes: the third q — 7 crossing (ki = 2, 
=9). The (adiabatic) states are labeled by letters and numbers indicate the (diabatic) resonance quantization k. Bigger dots 
mark on the quasienergy curves the points corresponding to the Husimi functions displayed. Lighter shades of gray between 
the level curves mean higher values of the Husimi functions. The levels of the curves is the same for all plots; the highest peaks 
in some of them are, as a result, out of range and appear as white areas. 
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FIG. 16: Detail ol Fig. 1 showing one ol the avoided crossings marked in Fig. 1 and the Husimi functions of the (eigen-)states 
undergoing this same avoided crossing: it is a wide crossing where the lower "adiabatic" state is made up of two states (B 
and C): the third q = 6 crossing (fci —2, k% — 8). The (adiabatic) states are labeled by letters and numbers indicate the 
(diabatic) resonance quantization k. Bigger dots mark on the quasienergy curves the points corresponding to the Husimi 
functions displayed. Lighter shades of gray between the level curves mean higher values of the Husimi functions. The levels of 
the curves is the same for all plots; the highest peaks in some of them are, as a result, out of range and appear as white areas. 
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FIG. 17: Detail of Fig. 1 showing one of the avoided crossings marked in Fig. 1 and the Husimi functions of the (eigen-)states 
undergoing this same avoided crossing: it is a wide crossing where the upper "adiabatic" state is made up of two states (A 
and C): the first q = 5 crossing (fei = 0, fa = 5). The (adiabatic) states are labeled by letters and numbers indicate the 
(diabatic) resonance quantization k. Bigger dots mark on the quasienergy curves the points corresponding to the Husimi 
functions displayed. Lighter shades of gray between the level curves mean higher values of the Husimi functions. The levels of 
the curves is the same for all plots; the highest peaks in some of them are, as a result, out of range and appear as white areas. 
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FIG. 18: Detail ol Fig. 1 showing one ol the avoided crossings marked in Fig. 1 and the Husimi functions of the (eigen-)states 
undergoing this same avoided crossing: it is a wwide crossing where a crossing state (B) maintains the same (diabatic) resonance 
quantum number k2 = 4: the first q — 4 crossing [kj = 0, k2 = 4). The (adiabatic) states are labeled by letters and numbers 
indicate the (diabatic) resonance quantization k. Bigger dots mark on the quasienergy curves the points corresponding to the 
Husimi functions displayed. Lighter shades of gray between the level curves mean higher values of the Husimi functions. The 
levels of the curves is the same for all plots; the highest peaks in some of them are, as a result, out of range and appear as 
white areas. 
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FIG. 19: The quasienergies of the states most significantly populated at every period of the microwave pulse for the case of 
Fig. 3 (F™ ax = 0.025). The (linear) size of each dot is proportional to the population on the corresponding level at that point, 
the minimum size corresponding to a 0.5% population. The same symbols as in Fig. 1 mark the avoided crossings related to 
classical secondary resonance island chains. 
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FIG. 20: The quasienergies of the states most significantly populated at every period of the microwave pulse for the case of 
Fig. 4 (_Fo" aa: = 0.029). The (linear) size of each dot is proportional to the population on the corresponding level at that point, 
the minimum size corresponding to a 0.5% population. The same symbols as in Fig. 1 mark the avoided crossings related to 
classical secondary resonance island chains. 
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FIG. 21: The quasienergies of the states most significantly populated at every period of the microwave pulse for the case of 
Fig. 5 (i?™ aa: = 0.041). The (linear) size of each dot is proportional to the population on the corresponding level at that point, 
the minimum size corresponding to a 0.5% population. The same symbols as in Fig. 1 mark the avoided crossings related to 
classical secondary resonance island chains. 
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FIG. 22: The quasienergies of the states most significantly populated at every period of the microwave pulse for the case of 
Fig. 6 (J?™ ax = 0.061). The (linear) size of each dot is proportional to the population on the corresponding level at that point, 
the minimum size corresponding to a 0.5% population. The same symbols as in Fig. 1 mark the avoided crossings related to 
classical secondary resonance island chains. The crosses mark crossings with p = 2 and -from top to bottom- q = 9, 9 and 10. 
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FIG. 23: Comparison of the k — — 9 curves obtained from eq. (B9) using the off diagonal matrix elements (BIO) (small dots) 
with the classical prediction (big dots) W.K.B. quantized according to the free atom quantum number n outside the separatrix 
and according to the resonance quantum number k inside. Both quantizations are shown at the separatrix itself, connected by 
a segment of straight line. (From Ref. [16]) 



